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ABSTRACT 



A digital filter representation of the ASQ-81 
magnetometer is derived from the s-plane transfer functions 
of the system through the use of a bilinear transformation. 
A FORTRAN computer program is written which applies this 
representation to time-sampled total magnetic field data in 
order to obtain a time series representation of ASQ-81 
filtered total field. A series of simulations and a field 
experiment are conducted which verify the program output. 
Applications of this program include usage in conjunction 
with geomagnetic field data in order to produce a new data 
set representative of geomagnetic noise observed by Navy MAD 
(Magnetic Anomaly Detection) aircraft with the potential to 
investigate techniques of reducing geomagnetic noise in MAD 
aircraft . 
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I. INTRODUCTION 



The detection and location of submarines (and other 
magnetic bodies) through the discrimination of changes or 
anomalies in the Earth's magnetic field is called Magnetic 
Anomaly Detection or MAD. In this technique, a magnetometer 
measures the magnitude of the Earth's magnetic field and 
provides an indication of that magnitude, or, more usually, 
an indication of changes in the magnitude of the Earth's 
field. These changes, or anomalies, can indicate the 
presence of magnetized bodies which may or may not be a 
submarine . 

The magnetometer currently in use in the United States 
Navy for use in this MAD process is the AN/ASQ-81 
metastable helium vapor total field magnetometer. 

Research is currently being conducted at the Naval 
Postgraduate School in Monterey, California, in various 
aspects of the applications of magnetometers, including 
Magnetic Anomaly Detection (MAD). Within the context of 
this research, magnetic field measurements are made through 
the use of sets of wire wound coils vice any specific 
magnetometer or magnetic detecting system. The data 
collected through the use of these coils is evaluated and 
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processed in a variety of methods for different project 
goals . 

This thesis project is designed to produce an acceptable 
alternative to the physical presence of an experimental 
AN/ASQ-81 magnetometer at the postgraduate school by 
allowing the determination, in conjunction with other 
research in progress, of the output of the AN/ASQ-81 
magnetometer from the data collected from the school's 
measurement coils. It is hoped that this will assist future 
research projects as, for example, in allowing a 
determination of environmental noise of such characteristics 
as to affect the AN/ASQ-81 magnetometer operationally with 
the eventual goal of providing an environmental noise index 
or a system of removing such noise from the magnetometer- 
detection system. 
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II. GEOMAGNETICS REVIEW 



A. EARTH'S MAGNETIC FIELD 

1 . Constituents of the Geoma gnet ic Field 

The most common method of specifying the consituting 
parts of the geomagnetic field is to divide the field in 
terms of distance from the center of the Earth. This 
method results in three classifications: internal . 
crustal . and external . [Ref. 1] 

The internal field originates in the core region and 
is the most stable field, containing only extremely low 
frequency temporal variations. The crustal, or anomalous, 
field arises from modifications made on the internal field 
by materials and structures in the Earth's crust. These 
variations are not constant with regard to spatial 
locations, and comprise part of what is known as geological 
variations. The external field is the most dynamic and 
arises from many sources, including the interaction between 
the solar wind and the Earth's magnetic field. 

In addition to this method of defining the Earth's 
magnetic field is the method of time variations. This 
method consists of considering that part of the field which 
varies with periodicities greater than about one year as the 
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steady field and everything else as the variation field . 



[Ref. 2] 

The steady field consists of the internal field , also 
refered to as the main field . Slow variations of the main 
field with periods of years or longer are referred to as 
secular variations . 

There are various elements that contribute to the 
geomagnetic field, some of which are external to the Earth's 
surface. External contributions make up only a small part 
of the steady field, but play an important role in the 
variation field. These external sources include current 
systems in the Earth's upper atmosphere affected by solar 
electromagnetic radiation and gravitation, solar corpuscular 
radiation and the interaction of solar plasma with the main 
field, and the effect of the solar inteplanetary field. 
[Ref. 3] 

The geomagnetic field changes with time. As 
previously mentioned, very slow variations in the main field 
with periods of on the order of years to thousands of years 
are referred to as secular variations. Secular variations 
are caused by a variation in the strength or orientation of 
the Earth's center dipole. 

Other time variations of the field can be categorized 
into quiet variation fields and disturbed variation fields . 
Disturbed variation fields include geomagnetic 
micropulsations , which are of particular interest to 
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operational forces as these can mask target signatures and 
are therefore a source of noise to MAD sensors. 

Quiet variation fields are those which are not due to 
disturbances in the interplanetary environment and which 
vary slowly and regularly. [Ref. 3] 

Disturbed variation fields are geomagnetic field 
variations that appear to be the result of interplanetary 
environmental changes and do not posess a simple 

periodicity. These variations include ionospheric 

disturbances, the aurora, geomagnetic storms, and 
geomagnetic micropulsations. 

2 . Elements of the Magnetic Field Vector 

The geomagnetic field vector is characterized at any 
point by its direction and magnitude. This is commonly 
accomplished through a system of coordinates as shown in 
Figure 2.1. The field is measured in terms of local 
coordinates with respect to true North. [Ref. 3] 

The various coordinates are referred to as ma gnetic 
elements and are defined as follows: 

Total field intensity (the symbol F^ is 

sometimes also used, as in this figure.) 

]1: Horizontal component 

Northward, or NorthSouth component 
X: Eastward, or EastWest component 

Downward, or Vertical component 
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Figure 2.1 Magnetic Field Elements [Ref. 4], 



Declination or magnetic variation 

This is the angle between IL 

is measured positive eastward. 

^ Inclination or Dip Angle. 

This is the angle between ]1 and ^ (or 
and is measured positive downward. 
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III. THE AN/ASQ-81 MAGNETOMETER 



A. DESCRIPTION OF SYSTEM OPERATION 

The Magnetic Anomaly Detecting set currently in use in 
the U S Navy is the AN/ASQ-81 magnetometer. This set is 
used to locate and classify submerged submarines by sensing 
disturbances in the Earth's magnetic field (anomalies) 
caused by the presence of the magnetic mass of the 
submarine. The disturbance of the Earth's field is detected 
by the magnetometer, processed through filtering circuits, 
and amplified. The output signal of the magnetometer is 
displayed on a chart recorder for interpretation by an 
operator . 

The magnetic detecting set is a metastable helium vapor 
magnetometer. The operation of the magnetometer is based on 
the light absorbtion properties of helium gas subjected to 
certain light stimulus (optical pumping), radio frequency 
excitation, and the Earth's magnetic field. The 
magnetometer consists of a helium lamp, lens and polarizer 
to generate a beam of polarized light radiation. This 
focused and polarized light beam is directed through a 
helium absorbtion cell to an infrared (IR) detector. Some 
of the helium gas in the absorbtion cell is maintained in a 
metastable state by application of VHF excitation. 
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Earfh*! mo^n«tic fi«ld vector 
/ 




/ 



Figure 3.1 : Metastable Helium Magnetometer [Ref«5]. 

The Earth’s magnetic field imposes a magnetic force upon 
the excited helium vapor atoms to force the atoms into one 
of three energy sublevels. This is called the Zeeman 
effect. The rate or frequency of atomic precession caused 
by this effect is called the Larmor frequency. A helium 
lamp is used to optically pump the atoms in the absorbtion 
cell, with the result that the polarized light energy 
passing through the absorbtion cell will polarize 
(magnetize) the helium atoms in the absorbtion cells by 
selectively pumping the Zeeman levels of the energy of the 
helium atoms in the cell. The magnetization direction is 
determined by the polarization of the photons from the 
helium lamp . 
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RF energy is then introduced to the absorbtion cell in 
the form of an additional magnetic field imposed through the 
use of coils oriented perpendicular to the precessed 
polarized helium atoms in the absorbtion cell and energized 
by a variable frequency RF oscillator. The RF oscillator is 
tuned to the Larmor frequency, which results in 
depolarization of the atoms. The atoms attempt to equally 
repopulate the Zeeman energy levels. However, the helium 
lamp is still beaming polarized light energy into the 
absorbtion cell, causing the atoms to absorb light energy 
and rise to an excited energy levbel. This absorbtion of 
light energy is detected through the use of an infrared 
detector. The RF oscillator frequency producing maximum 
light absorbtion is called the resonant frequency, and is 
determined through the use of a servo loop from the infrared 
detector to the RF variable frequency oscillator. 

Therefore, any change in the Earth’s magnetic field 
intensity will result in a change in the Larmor frequency of 
the helium atoms in the helium absorbtion cell. This new 
Larmor frequency will be detected by the ASQ-81 
magnetometer. Since the gyromagnetic ratio of helium is 
28.024 HZ per gamma, this detection of the resonant 
frequency provides a measurement of the Earth’s magnetic 
field intensity at any given time. A change in the Earth’s 
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magnetic field intensity could signal the presence of a 
submerged submarine. 

The output resonant frequency developed by the 
magnetometer is converted to a proportional output voltage 
which is filtered through the Magnetic Anomaly Detection 
(MAD) bandpass filters for environmental noise reduction 
and utilized to drive a chart recorder for observation by 
an operator. [Ref. 6] 

B. TRANSFER FUNCTIONS 

Transfer functions for the AN/ASQ-81 filters were 
obtained from the manufacturer of the AN/ASQ-81 detecting 
set, Texas Instruments of Dallas, Texas. These transfer 
functions are listed in Appendix A and are in the form of 
H(s), that is, the frequency domain, or S domain, where 
S = jw. The s-domain representation for transfer functions 
is routinely utilized to express output system 
characteristics for given system inputs. As the S domain 
representation is not utilzed further in this discussion, it 
will not be further explained. 

As the output signal of the ASQ-81 magnetometer is 
filtered through a fixed high-pass system, then through a 
selectable low pass system and a selectable high pass system 
(as shown in Figure 3.2 below), the transfer funtions are 
listed in this order. 
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Figure 3.2 ; Signal Flow Diagram for ASQ-81 

Further discussion will be made of the selectable filters 
later . 
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IV. DIGITAL FILTERING MODELLING OF SYSTEMS 



A. SEQUENCE REPRESENTATION OF TIME FUNCTIONS 
1 ♦ Signal Representation 

A signal can be defined as a function that conveys 
information, generally about the state or behavior of a 
physical system. Although signals can be represented in 
many ways, the information conveyed by the signal is 
contained in a pattern of variations of some form. Signals 
are represented mathematically as functions of one or more 
independent variables, one of the most common of which is 
time . 

The independent variable of the mathematical 
representation of a signal may be continuous or discrete. 
Continuous time signals are signals that are defined over 
continualy values of time and are therefore represented by 
continuous-variabled functions. Discrete time signals are 
defined at discrete time intervals and are therefore 
represented by functions whose independent variable(s) take 
on discrete values only. Discrete-time signals are 
represented as sequences of numbers . [Ref. 7] 

In addition to the fact that the independent 
variables can be either continuous or discrete, the signal 
amplitude can be either continuous or discrete. Digital 
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signals are those for which both time and amplitude are 



discrete. Analog signals are those for which both time and 
amplitude are continuous. 

Digital signal processing deals with transformations 
of signals that are discrete in both time and amplitude, 
usually represented by sequences of numbers. The nth number 
in the sequence 2L being processed is usually represented as 
X ( n ) . and is formally written as; 

x=[x(n)], -00< n < +O0 

In general, an arbitrary sequence can be expressed as 



x(n) 




x(k) d(n-k) 



where d(n-k) is the unit sampl e at time k. In other words, 
an arbitrary sequence may be expressed as a sum of scaled, 
shifted unit samples, where the scaling factor is equal to 
the amplitude of the sequence at that time. 

2 . Linear Shift-Invariant Systems 

A system is defined mathematically as a unique 
transformation or operator that maps an input sequence 
[x(n)] into an output sequence [y(n)]. This is denoted 
as : 

y(n) = T[x(n)] 
and is often depicted as in Figure 4.1. 
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I 

I 



■>“ y(n) 



! 

I 



x(n) 




Figure 4.1: Representation of Transformation of an Input 

Sequence to an Output Sequence. [Ref. 7] 

Classes of discrete time systems are defined by 

placing constraints on the transformation T[]. 

The class of linear systems is defined by the 

principle of superposition. If y (n) and y (n) are the 

1 2 

responses when x (n) and x (n) are the inputs, then a system 

1 2 

is linear if 

T[ax (n) + bx (n)] = aT[x (n)] + bT[x (n)] 

12 12 

= ay (n) + by (n) 

1 2 

for any arbitrary constants a and b. This, together with 

the concept of representing a sequence by a sum of delayed 

and scaled unit-sample sequences, suggests that a linear 

system can be characterized by its unit-sample response. 

Specifically, let h (n) be the response of the system to 

k 

d(n-k), a unit sample occurring at n=k . Then 

OO 



y(n)=T[ 7 x(k) d(n-k)] or, 

L I 

k= - OO 
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oo 



oo 



y(n) 




x(k) T[d(n-k)] 




x(k) 



k= -oo 



h (n) 
k 



Thus the system response can be expressed in terms of the 
response of the system to d(n-k). 

The class of shift invariant sy st ems is characterized 
by the property that if y(n) is the response to x(n), then 
y(n-k) is the response to x(n-k), where k is a positive or 
negative integer. When the index n is associated with time, 
shift-invariance corresponds to time-invariance . The 
property of shift invariance implies that if h(n) is the 
response to d(n), then the response to d(n-k) is simply 
h(n-k). Therefore 

y(n) = 

k= -OO 



X ( k ) h( n-k) 



and any linear shift-invariant system is completely 
characterized by its unit-sample response h(n). 

A subclass of linear shift-invariant systems consists 
of those systems for which the input x(n) and the output 
y(n) satisfy an Nth-order linear constant-coefficient 
difference equation of the form 



I 



a y(n-k) 
k 



Z 



b X ( n-r ) 
r 
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If the assumption of causality is made about the 
system, a linear difference equation provides an explicit 
relationship between the input to the system and the output 
of the system. This can be seen by rewriting the previous 
equation as 



y(n) 



Z 



M 



c y(n-k) + 
k 




d x(n-r) 
r 



where c =-a / a and c =b /a . 

k k 0 r r 0 

Thus the nth value of the output can be computed from 
the nth value of the input and the N and M past values of 
the output and input, respectively. The difference equation 
not only represents the system for theoretical purposes, 
but it may also serve as a computational realization of the 
system. The z-Transform makes use of this property to 
realize systems. 



B. THE z-TRANSFORM 

1 . Description of the z-Transform 

The z-tr ansf orm plays an important role in the 
analysis and representation of discrete-time linear shift- 
invariant systems. The z-transform, X(z), of a sequence 
x(n) is defined as 

X(z) = 

n= - CO 
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where z is a complex variable. This representation of the 
z-transform is referred to as the two-sided ^ transform . The 
one sided z-transform consists of the same summation for 
terms of n greater than or equal to zero. For the case 
that x(n)=0 for n<0, the one sided and two sided z 
transforms are equivalent. 

By expressing the complex variable z in polar form as 
z = re , the z-transform can be interpreted as the Fourier 
transform of x(n) multiplied by an exponential sequence. 
For r = 1, that is, for IzI = 1, the z-transform is equal to 
the Fourier tranform of the sequence. 

2 . The Bilinear Transformation 

The transfer functions of analog systems are most 
often expressed in terms of s = jw (see section III B.). 
This corresponds to the analog frequency response of the 
system. This analog frequency response can be "mapped”, 
that is, tranformed to the z-plane from the s-plane through 
the use of the bilinear transformation . The effect of 
utilizing the bilinear transformation is to convert a system 
transfer function in terms of the variable S into the system 
transfer function in terms of the variable z . The 
transformation itself is: 



and 



s 



-1 

2 1 - z 



-1 

T 1 + z 
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(2/T) + s 

z = 

(2/T) - s 

where T is the sampling period, that is, the time between 
data samples. 

Thus a transform can be made from one plane to the 
other. In this manner, the transfer function, H(z), of a 
system may be obtained. 

The bilinear tranf ormat ion equations may be shown to 
hold in general, and the use of this transformation may be 
shown to yield stable digital filters from stable analog 
filters [Ref. 7]. The bilinear transformation maps the 
imaginary jw axis in the s-plane onto a unit circle (of the 
region of convergence) in the z-plane, with the left half 
s-plane mapped onto the region inside the circle and the 
right hand (region of instability) s-plane mapped onto the 
region outside this circle [Ref. 8]. A complete discussion 
of the z-transform is available in several texts, some of 
which are listed in the Bibliography. 

C. THE DIGITAL COMPUTATIONAL ALGORITHM 

In implementing a digital filter on a digital computer 
such as the IBM 3033, the input-output relationship of the 
signals through the system being synthesized must be 

converted to a computational algorithm. The algorithm is 
specified in terms of a set of basic computations of 
elements. For the implementation of discrete- time systems 
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described by linear constant coefficient difference 
equations, such as the AN/ASQ-81, it is convenient to choose 
as these elements the basic operations of addition, delay, 
and multiplication by a constant. The computational 



algorithm 


for 


implementing the 


filter 


is then defined by 


a 


structure 


or 


network consisting of an 


interconnection 


of 


these basic 


operations. For a 


system 


transfer function 


of 


the form 




M 










k=0 Y(z) 

N X(z) 




k=l 



the difference equation relating input and output is easily 
written down directly from the system function and is given 
by N M 



y(n) 




a y(n - k) + 
k 




b x(n - k) 
k 



[Ref. 7] 



k=l 



k=0 



This difference equation can be interpreted directly as a 

computational algorithm in which the delayed values of the 

input are multiplied by the coefficients b , the delayed 

k 

values of the output are multiplied by the coefficients 

a , and the resulting products are added. It is now easy 
k 

to see the process to be followed in obtaining the 
computational algorithm for the AN/ASQ-81 magnetometer 
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transfer function. The z-transform of the system transfer 
function is obtained through the use of the bilinear 
transformation, and is then converted into a difference 
equation relating input and output signals, thence to a 
FORTRAN computer program. A table of z-transforms of system 
functions is included in Appendix B. 

In the FORTRAN computer program realization of the total 



system computational 


algorithm , 


each filter block 


is 


transformed 


into a 


separate 


difference equation 


and 


algorithm . 


This was 


done to enable a "building block" 


type 


approach to 


the program, and to 


minimize computational 


and 



roundoff errors. 

D. THE CASCADE FORM OF THE COMPUTATIONAL ALGORITHM 

Even though the direct form realization of the digital 
filter design may be perfectly satisfactory in a theoretical 
sense, it may be less than desirable in the context of 
realization through the use of a general purpose computer of 
fixed register length. The parameters of a digital filter 
are usually obtained with a high degree of accuracy, which 
results in a faithful realization of the desired system. 
When these parameters are quantized, as in a finite memory 
register within a computer, the frequency response of the 
resulting digital filter may differ appreciably from the 
original design. In fact, the quantized filter may fail to 



28 



meet design specifications although the unquantized filter 
does. [Ref. 7 \ 



The 


sensitivity 


of 


the filter 


response to 


errors 


in 


the 


filter 


parameters 


is 


dependent 


upon the structure 


of 


the 


filter 


realization . 




Therefore 


, in the 


event 


of 


an 



unacceptable change in the frequency response of the filter 
due to quantization errors, it is often possible to minimize 
the effect of these errors through an alternate filter 
realization structure. An alternate structure to the 
previoulsy discussed direct form realization is the cascade 
form realization. 

The direct form network structures were obtained directly 

from the system function H(z) written in the form of a ratio 

of sums. If this ratio is factored into a product of 

polynomials of the form 

[(N+l)/2] _2 

I , I 1 + B z + B z 

^ Ik 2k 



H(z) = A 



-1 -2 



k=l 



1 - a z - a z 
Ik 2k 



this product represents a general distribution of poles and 
zeros and suggests a set of structures consisting of a 
cascade of first and second-order subsystems. There is 
considerable freedom in the design of the subsystems, but it 
is best to realize the systems using a minimum of storage. 

The expression of H(z) in this form indicates the 
presence of poles and zeros in pairs. If poles and zeros 
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are not present in pairs, one of the coefficients B or a 

2k 2k 

will be zero as appropriate. An implementation of such a 
cascade structure with the use of minimum memory can be 
obtained through a direct form II realization of each second 
order subsystem using techniques similar to the direct form 
implementation utilized previously. A cascade realization 
of a sixth-order system, such as the ASQ-81 system, using a 
direct form II realization of each second order subsystem 
would appear as in Figure 4.2 below. 
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Figure 4.2: Cascade Structure With a Direct Form II 

Realization of Each Second Order Subsystem. 
[Ref. 7] 



There is, theoretically, considerable flexibility in the 
manner in which the poles and zeros are paired together and 
in the order in which the resulting second-order subsystems 
are cascaded. However, although all such pairings and 
orderings are equivalent for infinite-precision arithmetic, 
they may differ considerably in practice owing to finite 
word length effects of roundoff and truncation. 
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Figure 5.1: Output of First Stage Filter of Digital Filter 
Computer Program With Sinusoidal Input in Simulation. 

Unfortunately, the second stage output of the filter 
showed an instability within the program design, indicated 
by the output of the filter being a sinusoid of increasing 
magnitude, as indicated in Figure 5.2 below. 
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Figure 5.2: Output of Second Stage Filter Design With Input 

of a Sinusoid. 
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The stability of the third stage of the filter design 
was investigated by inputting the sinusoid directly to the 
third filter, and found to be stable. A check of the 
derivation of the equations, coefficients, and programming 
steps of the second (unstable) filter of the design failed 
to indicate the cause of the instability. 

Computation of the poles of the z transfer function, 
H(z), of the second stage of the filter confirmed the 
instability of the design. The poles were computed to be: 
0.92+0.12181, 1.07+0.13401, 0.8611, and 1.1564. Of 

these six poles, three lie outside the region of convergence 
for the z-plane, that is, within the unit circle discussed 
previously in Chapter IV. 

The second stage of the filter was therefore 
redesigned using the cascade form of the direct form 
realization (direct form II), and tested in simulation. A 
copy of the software used in the simulation is enclosed in 
Appendix F. 

The output of all three filter stages of the program 
were stable, as indicated in Figures 5.3 through 5.7 below. 
The amplitude decrease and phase shift expected were 
observed. The "damped overshoot" of the second stage output 
is due to the fact that, for values of the input function 
prior to time zero in the simulation, utilized in the input- 
output signal difference equations for the filter, the input 
signal was set at 0. This resulted in an instantaneous 
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change of the input signal from 0 to the finite value 
introduced in the simulation at time 0+ . The ’’overshoot’’ of 
the* filter is the filter’s attempt to ’’match” this 
instantaneous jump in magnitude of the input signal. When 
the input signal to the filter in the simulation is zero 
at time zero, this overshoot effect does not occur. 
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Figure 5.4: Output of First Stage of Digital Filter Program. 
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Figure 5.5; Output of Second Stage of Digital Filter 

Program. 



35 
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Figure 5.6: Output of Third Stage of Digital Filter Program. 

This simulation was run with inputs of sinusoids of 
various frequencies in order to check the stability of the 
filter design at frequencies throughout the operating range 
of the AN/ASQ-81 magnetometer. In all cases, the design was 
stable, and the expected amplitude changes and phase shifts 
occured . 

2 . Noiselike Inputs 

The simulation was also run with inputs consisting of 
a sinusoid of a frequency which should be passed through the 
AN/ASQ-81 added to sinusoids of frequencies which should 
have been filtered by the magnetometer and random noise. The 
filter performed as expected, with the sinusoid of a 
passable frequency passed by the filter, and spurious noise 
and sinusoids attenuated severely. The results of a 
simulation consisting of a sinusoid of passable frequency, a 
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filterable sinusoid, and uniformly distributed random noise, 
all of amplitude ^ 1, are presented in Figures 5.7 through 
5.10 below . 
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Figure 5.7: Input to Filter - 0.1 HZ Sinusoid, 10 HZ 

Sinusoid, Uniformly Distributed Random Noise of Amplitude 

~l" 1 • 




Figure 5.8: Output of First Filter Stage. 
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Figure 5.9: Output of Second Filter Stage. 




As can be seen, the digital filter program succeeds 
in filtering out random noise and signals of frequency 
components above the band pass of the magnetometer. 

In order to ensure that the digital filter 
representation of the magnetometer has the same amplitude 
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versus frequency characteristics of the AN/ASQ-81 
magnetometer, a simulation program was written which inputs 
sinusoids of varying frequencies and computes the Root Mean 
Square (RMS) value of the filter output and the signal 
input, then computes the decibel (dB) attenuation of the 
filter at that frequency. A copy of this program is 
included in Appendix G. A plot was made of the dB 
attenuation versus frequency for the filter and compared 




PREQUENCY - HZ 



Figure 5.11; Plot of Attenuation Versus Frequency for 
Sinusoidal Inputs for Digital Filter and Analog Filter. 
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with the measured frequency performance of the AN/ASQ-81 
magnetometer, which was supplied by Texas Instruments, Inc., 
and does not include the effects of the fixed high pass 
filter. Consequently, the data shown in Figure 5.11 is a 
comparison of the data supplied by Texas Instruments and the 
output of the test program, which also does not include the 
fixed high pass filter. As can be seen, the performance of 
the filter is extremely similar to that of the magnetometer 
itself . 

3 . Anderson Function Simulations 

The next step in the simulation phase was the 
introduction to the filter of Anderson function simulations. 
The shape of the signal amplitude of the output of a 
magnetometer passing through the sphere of influence of a 
magnetic anomaly (submarine) is a function of the dip angle 
of the geomagnetic field, the magnetic heading of the track 
of the magnetometer (or the aircraft), the magnetic 
heading of the anomaly (submarine) dipole, and the 
lateral range between the magnetometer (aircraft) and the 
anomaly. Anderson functions [Ref. 9] are mathematical 
representations of three basic components of signals which, 
when taken in various linear combinations, describe the 
shape of these anomaly signals. The equations for the 
Anderson functions are: 
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2 5/2 

(First Anderson Function) f = 1/(1 + B ) 

0 

where B = (velocity)x( time ) 

i 

range at CPA 

or, a dimensionless parameter- defined as the distance 
traveled along the magnetometer (aircraft) track divided by 
the slant range at closestpoint of approach (CPA) 

(Second Anderson Function) f = B x f 

1 0 2 

(Third Anderson Function) f =Bxf =B xf 

2 1 0 

The Anderson functions were introduced into the 
filter program in a noise-free signal environment in order 
to observe the output signal and ensure that it was a "MAD- 
like" signal. A rigorous determination of the actual output 
signal would have been extremely difficult to obtain, so a 
comparison was made with the output of a computer simulation 
program provided to NPS by Mr. Joe Rice of Texas 
Instruments. When the sampling rate of the program was 
adjusted to equal that of the Texas Instruments program, 8 
HZ, the two program outputs were observed to be very 
similar. The Anderson function simulation inputs and 
outputs of the program are depicted in Figures 5.12 through 
5.18. The Texas Instruments program outputs were obtained 
in the form of time series plots of discontinuous data 
points, and were therefore not conducive to replotting for 
comparison . 
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Figure 5.12: First Anderson Function Input. CPA at Time 

2.5 Minutes . 




Figure 5.13: Filter Output for First Anderson Function 

Input . 
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Figure 5.14: Second Anderson Function Input. CPA at Time 

2.5 Minutes. 
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Figure 5.16: Third Anderson Function Input. CPA at Time 

2.5 Minutes. 
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Figure 5.17: Filter Output for Third Anderson Function. 

The filter output for all three Anderson function 
inputs did appear to be "MAD-like" signals, and did closely 
resemble the simulation output obtained from Texas 
Instruments, Inc. 

4 . Impulse Function Response 

The response of the filter program was also observed 
when the input was a unit impulse function. Again, the 
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Figure 5.18: Impulse Response of Filters. 

of the Texas Instruments response are due to the translation 
of the output plot supplied to this plot. The plot supplied 
by Texas Instruments was, again, discontinuous points of 
poor resolution, and it was necessary to interpolate values 



45 



in order to generate Figure 5.18. This resulted in the 
broken appearance of the plot. Even so, the similarity of 
the outputs can be observed. 

B. EQUIPMENT SETUP 

Following the simulation phase of the experiment, actual 
magnetic field measurements were introduced to the filter in 
order to test the response of the filter. Magnetic field 
measurements were made at the La Mesa field test site near 
the Naval Postgraduate School in Monterey. The output of an 
AN/ASQ-81 magnetometer, a Schonstedt magnetic field sensor, 
and the school's coil sensor, oriented along the Earth's 
magnetic field, were pulse code modulated (PCM) and 
transmitted via VHF radio to recording devices at the 
Postgraduate school. The recording of a two hour long data 
collection period was transferred to digital data tape for 
use by the school's IBM3033 general purpose mainframe 
computer . 

In the first test of the digital filter program, the 
output of the Schonstedt sensor, which represents 
fluctuations of the Earth's total field, was used as the 
input to the computer program. A comparison of the output 
of the computer program, with this approximation to the 
total field fluctuations as input, to the output of the 
AN/ASQ-81 should provide an indication of the proper 
functioning of the computer filter program. The results of 
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the test are shown in Figures 5.19 through 5.21 on the pages 
following. Figure 5.19, the Schonstedt sensor output, shows 
several instances of PCM dropouts, that is, occasions where 
the pulse code modulation signal was not correctly read by 
the computer for some reason. At such occurances, the data 
point value used by the computer is a random number and does 
not reflect the true value of the data. The problem with 
these PCM dropouts is that the computer does not recognize 
them as invalid data points and will use them in 
computations. This can (and does) cause problems in the 
computation of Fourier transforms, spectral characteristics, 
etc. Additionally, this will also impact the proper 
functioning of the digital filter program which is the 
subject of this thesis. PCM dropouts are visible at times 
6, 8, 142, and 220 through 226 seconds on the plot of the 
Schonstedt sensor output. An examination of Figure 5.20, 
the filter program output, reveals the programs attempt to 
"follow” these PCM dropouts. It should be recalled that 
previous simulations indicated the filter's tendency to 
"follow" sudden changes in the input signal, with a 
relaxation time required for the filter to steady out. This 
effect is apparent in the output of the filter program at 
times corresponding to those of the PCM dropouts in the 
Schonstedt sensor's time series plot. It can be seen that 
this overshoot tendency resulted in an output significantly 
different from the actual AN/ASQ-81 output at these times. 
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Figure 5.19: Schonstedt Coil Time Series Output 
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Figure 5.20: Program Time Series Output Witli Schonstedt Coil as Input 
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Figure 5.21: ASQ-81 Time Series Output 



If the PCM dropout induced differences are neglected, it can 
be seen that the shape of the output of the filter program 
is remarkably similar to that of the AN/ASQ-81 magnetometer, 
although noisier. Note that the output of the AN/ASQ-81 
magnetometer exceded the maximum voltage amplitude which the 
pre-amplifiers of the data collection system were able to 
handle and resulted in a truncated signal from time 40 to 
time 60 seconds. It can still be seen, however, that the 
filter program output is very similar to the time signal 
which would have been displayed without this truncation. 

It should be noted that the amplitudes of the time series 
signals of the program output and the AN/ASQ-81 magnetometer 
differ considerably. In the case of the program output, the 
peak amplitudes are on the order of 1.4 nanoteslas along the 
vertical scale, while the peak amplitudes of the output of 
the AN/ASQ-81 magnetometer are on the order of 0.7 
nanoteslas along the vertical scale. This is because the 
input signal to the filter program is an approximation to 
the total field difference time series signal, and some 
amplitude difference could reasonably be expected. The 
intent of this initial test was to investigate the output 
time series shape, and an exact correlation was not 
expected. It is worth noting that the digital filter 
program will perform its function on any time series signal, 
regardless of units. This means that a signal may be 
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operated upon either before or after conversion from 
whatever units it was originally measured to magnetic field 
strength units. 

Therefore it appears that the digital filter program is 
functioning properly. When a close approximation to the 
fluctuations of the total field time series signal is used 
as the input to the computer program, the output of the 
program is similar to the time series output of an AN/ASQ-81 
magnetometer . 

The final stage in the testing process was a conversion 
of the time series output voltage signal of the coil antenna 
sensor, which was aligned along the Earth's magnetic field, 
into a total field fluctuation time series representation 
for the same time period as before, and then to use this as 
the input to the digital filter program. A comparison of 
the resultant time series output of the program with the 
actual AN/ASQ-81 magnetometer output would validate the 
proper functioning of the program. 

Conversion of the time series antenna sensor output 
voltage signal into total field fluctuations in nanoteslas 
was accomplished through the use of a computer program 
designed by Capt. Kurt Stevens, USAF, a student at the Naval 
Postgraduate School, as his Master's thesis [Ref 10]. The 
output voltage time series is stored in an array, then a 
Fourier transform is performed on the stored data, resulting 
in the Fourier spectrum of the data. This spectrum is 
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corrected for the characteristics of the coil antenna sensor 



to obtain the Fourier spectrum of the total field data. A 
reverse Fourier transform gives the time series signal for 
total magnetic field in nanoteslas. 

This time series signal was used as the input to the 
digital filter program and compared with the output of the 
AN/ASQ-81 magnetometer. Figures 5.22 through 5.24 show the 
raw coil antenna data, the total field time series data, and 
the program output time series for a 6 minute period of the 
test. Figure 5.22 shows the raw coil antenna data series. 
The number of PCM dropouts should be noted, as these will 
influence the performance of the filter program. Figure 
5.23 shows the computed total field time series. Note that 
the PCM dropouts evident on the raw time series plot are 
evident on the computed total field time series plot also, 
and thus inputted to the filter program as valid data 
points. Additionally, there are two "jumps" in the plot of 
total field fluctuation (Figure 5.23) which are also 
inputted to the filter program as valid data points. These 
"jumps" are located at 128 and 256 seconds and are caused by 
the method of processing blocks of data for the conversion 
to total field fluctuation. A block of 128 seconds of data 
is processed at a time, and the results of each block are 
stored in an array. This results in a slight amplitude 
difference between the last data point of one block and the 
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Figure 5.22: Raw Coil Antenna Time Series Output 
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Figure 5.23; Coii Antenna Difference Field Time Series Output 
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Figure 5.24: Program Time Series OuLpuL With Coil Anteuna Difference Field 

Time Series as Input. 



first data point of the next block of data. This slight 
difference is manifested as a signal jump. 

A comparison of Figures 5.24 (program output) and 5.21 
(AN/ASQ-81 output) show that the filter program gives a time 
series output very similar to that of the actual 
magnetometer. The first 20 seconds of the program output is 
somewhat dissimilar to that of the AN/ASQ-81, due either to 
the initial ’’start up” delay of the filter program or to 
distortion of the total field fluctuation time series. 
There is a PCM dropout at time 11 seconds which contributed 
to the distortion. 

Following this, however, it can be seen that the program 
output is very similar to that of the magnetometer, except 
at 128 and 256 seconds, which show the effects of the false 
signal jumps caused by the total field fluctuation 
conversion. There is also a noticeable four to five second 
time delay between the AN/ASQ-81 output and that of the 
filter program. As this delay is not evident in a 
comparison of the AN/ASQ-81 output and that of the filter 
program with the Schonstedt sensor as the input, it can be 
inferred that this time delay is caused either by the 
program which converts the raw coil data to total field 
fluctuation data, or by a phase (and hence time) change of 
the voltage signal due to the coil sensor itself. A 
comparison of the raw coil data in Figure 5.22 to the 
converted coil data in Figure 5.23 indicates no time shift. 
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and hence the deduction can be made that there is a time 



delay inherent within the coil sensor itself. 

Other than the differences of the four to five 
time delay and the distortions caused by the false 
jumps, the program output is extremely similar to the 
of the AN/ASQ-81 magnetometer. 



second 

signal 

output 
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VI, CONCLUSIONS 



The intent of this thesis was to design and test a 
digital filter computer program which would, when given a 
time series input of fluctuations in the total magnetic 
field, deliver an output time series representation of the 
output of an AN/ASQ-81 magnetometer. This purpose has been 
realized . 

The computer program contained in Appendix I has been 
proven to output a time series signal which is very similar 
to that of the magnetometer. The major limitations of the 
output signal are a finite time delay of about five seconds 
between the AN/ASQ-81 magnetometer signal and the output 
signal of the program, a sensitivity of the program to false 
data points such as those caused by PCM dropouts and false 
signal jumps caused by processing large data blocks, and the 
inherent limitations of the program caused by its dependence 
on the use of digital data tapes and the IBM 3033 mainframe 
computer . 

The five second time delay is not considered to be an 
important limitation to the program, as it was intended as a 
research tool for programs currently in progress at the 
Naval Postgraduate School. Instances where this time delay 
might become important would be in areas of simultaneous 
comparison of the program output signal with an actual 
magnetometer, in target location algorithms using time 
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delays, or in correlation studies between different sensors. 
In correlation studies using coil sensors, the effects of 
the time delay would cancel out, as all coil outputs would 
be similarly delayed. In target location algorithms, the 
target location errors due to the time delay could be 
adjusted for simply while in computer simulation, and flight 
testing could not reasonably be accomplished without the use 
of an actual magnetometer as the sensor. Lastly, in a 
comparison of the program output with an actual sensor, the 
time delay can, again, be compensated for. In short, these 
limitations are not considered excessive, especially as the 
apparent cause for the delay is not the filter program. 

In the primarily intended purpose of the filter program, 
magnetic noise studies, the time delay is not considered to 
be a problem. 

The problem of false data points caused by PCM dropouts 
and signal jumps due to conversion to total field 
fluctuations is more serious. False data points cause 
inaccuracies in the output time series and could adversely 
affect later projects. Unfortunately the PCM dropout 
problem is one which is endemic to the data collection 
system presently being used at the postgraduate school, and 
not to the filter program itself. It is imperative that 
users of this program are aware of the PCM dropout problem 
and of the effects it may entail upon their specific 



60 



research. A large number of PCM dropouts in a time series 
could render that series unusable. Similarly the case of 
the false data jumps caused by conversion to total field 
fluctuations is not within the filter program. Further 
investigation of this problem is necessary in order to 
eliminate it . 

The last problem, that of reliance upon the digital data 
tape/IBM 3033 computer system, is, like the PCM problem, one 
which is not endimic to the filter program but rather to the 
data collection system being used. A change of data 
collection system may, at some future time, remove the 
reliance upon the PCM/digital tape/IBM 3033 data system (and 
hence too the data block conversion problem which results in 
falso data jumps), but this is unlikely at this time. Users 
should be aware of this dependence and of possible effects 
upon specific research projects. 
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APPENDIX A 



AN/ASQ-81 FILTER TRANSFER FUNCTIONS 



Fixed High Pass Transfer Function ; 

2 

80 S 

H(S) = 



80 S + 20 S + 1 

Selectable High Pass Transfer Functions: 

A. 0.04 ^ H(S) = 

2 

40.82834 S 



X 

40.82834 S + 12.52096 S + 1 45.28317 



57.576688 S 



57.576688 S^ + 7.41498 S + 1 



2 

45.28317 S 
s"^ + 11.00999 S + 1 



B 0.06 HZ H(S) = 

2 2 
18. 14591 S 20.12587 S 



18.14591 S^ + 8.34727 S + 1 20.12587 S" + 7.33999 S + 1 

2 

25.58964 S 

X ^ 

25.58964 S“ + 4.94332 S + 1 



C. 0.08 HZ H(S) = 

2 

10.20708 S 



10.20708 S + 6.26045 S + 1 

2 

14.39417 S 



14.39417 S + 3.70749 S + 1 



2 

11.32080 S 



11.32080 S^ + 5.50500 S + 1 
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D. 0.10 HZ 



H(S) 



2 

6.53253 S 



6.53253 S + 5.00836 S + 1 

2 

9.21227 S 



9.21227 + 2.96599 S + 1 



2 

7.24531 S 

7.24531 + 4.40400 S + 1 



Selectable Low Pass Transfer Functions 



A. HZ 


H(S) = 






1 


1 


0.3143 S ^ + 


1.0741 S + 1 


0.2501 + 0.6209 S + 1 


B. 0^ HZ 


H(S) = 






1 


1 


0.07858 + 


0.53706 S + 1 


0.06252 + 0.31044 S + 1 


C. 0.6 HZ 


H(S) = 






1 


1 

V' ^ .. .. . 



2 ^ 7 

0.03492 S + 0.35804 S + 1 0.02779 S“ + 0.20696 S + 1 
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APPENDIX B 



AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS FOR DIRECT 

FORM I REALIZATION 

For fixed high pass filter : 

-1 -2 
BFHPO + BFHP1*Z + BFHP2»Z 

H(Z) = -_2 

1 - AFHP1*Z - AFHP2»Z 

where BFHPO, BFHPl, BFHP2 , AFHPl, , AFHP2 are constants 

tabulated in Appendix D. 

For selectable high pass filter ; 

-1 -2 -3 

BSHPO + BSHP1*Z + BSHP2*Z + BSHP3*Z 

+BSHP4*Z~^ + BSHP5*Z"^ + BSHP6*Z 

H(Z) -P 

1 - ASHP1*Z - ASHP2*Z - ASHP3*Z 

- ASHP4*Z'^ - ASHP5*Z - ASHP6»Z~^ 

where, for low frequency cutoff of 0 . 04 HZ : 

BSHPO = 0.99471378 



BSHPl 


= 


-5.9682827 


ASHPl 


= 


5.9894021 


BSHP2 


= 


14.920707 


ASHP2 


= 


-14.947051 


BSHP3 


= 


-19.894276 


ASHP3 


= 


19.894225 


BSHP4 


= 


14.920707 


ASHP4 


= 


-14.894327 


BSHP5 


= 


-5.9682817 


ASHP5 


= 


5.9472141 


BSHP6 


— 


0.99471372 


ASHP6 




-0.98945296 
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J"tL ^ 












For low 


frequency cutoff 


of 


0.06 


HZ; 






BSHPO = 


0.9920813 












BSHPl = 


-5.9524928 






ASHPl 


= 


5.9841070 


BSHP2 = 


14.881232 






ASHP2 


= 


-14.920650 


BSHP3 = 


-19.841643 






ASHP3 


= 


19.841528 


BSHP4 = 


14.881232 






ASHP4 


= 


-14.841757 


BSHP5 = 


-5.9524927 






ASHP5 


= 


5.9209919 


BSHP6 = 


0,99208212 






ASHP6 


= 


-0.98422128 


For low 


frequency cutoff 


of 


0.08 


HZ; 






BSHPO = 


0.98945806 












BSHPl = 


-5.9367483 






ASHPl 


= 


5.9788144 


BSHP2 = 


14.841871 






ASHP2 


= 


-14.894276 


BSHP3 = 


-19.789161 






ASHP3 


= 


19.788959 


BSHP4 = 


14.841871 






ASHP4 


= 


-14.789365 


BSHP5 = 


-5.9367476 






ASHP5 


= 


5.8948841 


BSHP6 = 


0.98945802 






ASHP6 


= 


-0.97901720 


For low 


frequency cutoff 


of 


0.1 


HZ; 






BSHPO = 


0.98684156 












BSHPl = 


-5.9210493 






ASHPl 


= 


5.9735244 


BSHP2 = 


14.802623 






ASHP2 


= 


-14.867941 


BSHP3 = 


-19.736831 






ASHP3 


= 


19.736516 


BSHP4 = 


14.802623 






ASHP4 


=r 


-14.737149 


BSHP5 = 


-5.9210491 






ASHP5 


= 


5.8688889 


BSHP6 = 


0.9868415 






ASHP6 


— 


-0.97384065 
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For selectable low pass filter ; 

-1 -2 
BSLPO + BSLP1*Z + BSLP2»Z 

H(Z) = --- 

1 - ASLP1*Z - ASLP2*Z ^ - 

where BSLPO, BSLPl, BSLP2, BSLP3, 

ASLP3, ASLP4 are constants tabulated 



-3 -4 

+ BSLP3*Z + BSLP4*Z 

ASLP3*Z"^ - ASLP4*Z 

BSLP4, ASLPl, ASLP2, 

in Appendix D. 
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APPENDIX C 



AN/ASQ-81 Z TRANSFORM FILTER TRANSFER FUNCTIONS DIRECT 

FORM II REALIZATION 



For fixed high 



H(Z) 

where BFHPO, 



pass filter ; 

-1 -2 
BFHPO + BFHP1*Z + BFHP2*Z 



1 - AFHP1*Z‘^ - AFHP2*Z“2 
BFHPl, BFHP2, AFHPl, AFHP2 are 



constants 



tabulated in Appendix D. 



For selectable high pass filter ; 

-1 -2 
1 _ 2*7 + 7 

H(Z) = ASHPl X -_-r X 

1 - ASHP3*Z ^ - ASHP2*Z 

-1 -2 -1 -2 
1 - 2*Z + Z 1 - 2*Z + Z 



1 - ASHP4*Z ^ - ASHP5*Z"^ 1 - ASHP6*Z ^ - ASHP7*Z 



where ASHPl, ASHP2, ASHP3, ASHP4, ASHP5, ASHP6, ASHP7 are 
constants and tabulated in Appendix D. 

For selectable low pass filter ; 

-1 -2 -3 -4 

BSLPO + BSLP1*Z + BSLP2*Z + BSLP3*Z + BSLP4*Z 

H(Z) -2 -5 

1 - ASLP1*Z - ASLP2*Z - ASLP3*Z - ASLP4*Z 

where BSLPO. BSLPl, BSLP2, BSLP3, BSLP4, ASLPl, ASLP2, 
ASLP3, ASLP4 are constants tabulated in Appendix D. 
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APPENDIX D 



Z TRANSFORM REALIZATION DIFFERENCE EQUATIONS 

With reference to Figures 3.2 and 4.2, the following 
difference equations are used to model the AN/ASQ-81 
magnetometer filter transfer functions. The input to the 
fixed high pass filter is called SIG(I), where I is the 
current data sample. The output of the fixed high pass 
filter, which is the input to the selectable high pass 
filter, is YO(I), and the output of the selectable high 
pass filter, the input to the selectable low pass filter, 
is called YPO(I). The output of the filter is called 
ASQ(I). (I-l) denotes a time delay of one sample, and so 
forth, and the symbol * denotes multiplication. 

For the fixed high pass filter ; 

YO(I)=BFHPO*SIG(I) + BFHP1*SIG(I-1) + BFHP2*SIG( 1-2 ) + 

AFHP1*Y0(I-1) + AFHP2*Y0(I-2) 

where : 

BFHPO = 0.9980499222938581 

BFHPl = -1.9960998445877161 AFHPl = 1.9960983216843922 

BFHP2 = 0.998049922238581 AFHP2 = -0.9961013674910398 

For the selectable high pass filters ; 

XI(I) = ASHP1*Y0(I) + ASHP2*XI(I-2) + ASHP3*XI ( 1-3 ) 

XII(I) = XI(I) + XI(I-2) - 2*XI(I-1) 

XIII(I) = XII(I) + ASHP4»XIII(I-1) + ASHP5*XIII(I-2) 

XIV(I) = XIII(I) - 2-»-XIII(I-l) + XIII(I-2) 



68 



XV(I) = XIV(I) + ASHP6*XV(I-1) + ASHP7*XV(I-2) 



YPO(I) ^ 


= XV(I) - 2*XV(I-1) + 


XV(I-2) 




For the 


low frequency cutoff 


at 0.04 


HZ: 


ASHPl = 


0.994713789347288 


ASHP2 = 


-0 


ASHP3 = 


1.9952137256322473 


ASHP4 = 


1 


ASHP5 = 


-0.9962082013015601 


ASHP6 = 


1 


ASHP7 = 


-0.9979897681491607 






For the 


low frequency cutoff 


at 0.06 


HZ: 


ASHPl = 


0.9920821277199399 


ASHP2 = 


-0 


ASHP3 = 


1.9928247245317958 


ASHP4 = 


1 


ASHP5 = 


-0.9943177045263504 


ASHP6 = 


1 


ASHP7 = 


-0.9969861717656565 






For the 


low frequency cutoff 


at 0.08 


HZ: 


ASHPl = 


0.9894580558875787 


ASHP2 = 


-0 


ASHP3 = 


1.9904384565912725 


ASHP4 = 


1 


ASHP5 = 


-0.9924307799269113 


ASHP6 = 


1 


ASHP7 = 


-0.9959835860201267 






For the 


low frequency cutoff 


at 0.10 


HZ: 


ASHPl = 


0.9869415560096681 


ASHP2 = 


-0 


ASHP3 = 


1.9880549117849344 


ASHP4 = 


1 


ASHP5 = 


-0.9905474442556225 


ASHP6 = 


1 


ASHP7 = 


-0.9949820174650785 






For the 


selectable low pass 


filters ; 




ASQ(I) = 


= ASLP1*ASQ(I-1) + ASLP2*ASQ(I 


-2) 



+ ASLP4*ASQ(I-4) + BSLPO*YPO(I) + 
+ BSLP2*YP0(I-2) + BSLP3*YP0(I-3) 



.9952196910157882 

.9962028201103847 

.9979855321466768 

.9928381306174365 

.9943056083414792 

.9969766455640039 

.9904622611337722 

.9924092959984783 

.9959666590811389 

.9880929619779491 

.9905139074397893 

.9949555824611562 

+ ASLP3*ASQ(I-3) 
BSLP1*YP0(I-1) 

+ BSLP4*YP0(I-4) 
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For the 
BSLPO = 
BSLPl = 
BSLP2 = 
BSLP3 = 
BSLP4 = 
For the 
BSLPO = 
BSLPl = 
BSLP2 = 
BSLP3 = 
BSLP4 = 
For the 
BSLPO = 
BSLPl = 
BSLP2 = 
BSLP3 = 
BSLP4 = 



high frequency cutoff 


at 0.2 


HZ: 


0.0000000452616229 






0.0000001810464917 


ASLPl = 


3.9082436339591027 


0.0000002715697375 


ASLP2 = 


-5.7285022156249328 


0.0000001810464917 


ASLP3 = 


3.7321935213310065 


0.0000000452616229 


ASLP4 = 


-0.9119356638511430 


high frequency cutoff 


at 0 . 4 


HZ: 


0.0000006918001209 






0.0000027672004837 


ASLPl = 


3.8173771378993420 


0.0000041508007256 


ASLP2 = 


-5.4670046844062743 


0.0000027672004837 


ASLP3 = 


3.4812554457127576 


0.0000006918001209 


ASLP4 = 


-0.8316389680077603 


high frequency cutoff 


at 0 . 6 


HZ: 


0.0000033463317975 






0.0000133853271900 


ASLPl = 


3.7274299052305002 


0.0000200779907850 


ASLP2 = 


-5.2152772583906819 


0.0000133853271900 


ASLP3 = 


3.2462216641520829 


0.0000033463317975 


ASLP4 = 


-0.7584278523006615 
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APPENDIX E 



DIGITAL SOFTWARE FOR SIMULATION - DIRECT FORM (SINUSOIDS AS INPUT) 
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APPENDIX F 

DIGITAL SOFTWARE FOR SIMULATION - CASCADE FORM (SINUSOIDS AS INPUT) 
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APPENDIX G 

DIGITAL SOFTWARE FOR COMPUTATION OF SYSTEM AMPLITUDE VERSUS FREQUENCY 
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APPENDIX H 

DIGITAL SOFTWARE FOR SIMULATION (ANDERSON FUNCTIONS AS INPUT) 
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